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Abstract 

We present a combined theoretical and experimental study of the buckling of a 
thin film wrapped around a sphere under the action of capillary forces. A rigid sphere 
is coated with a wetting liquid, and then wrapped by a thin film into an initially 
cylindrical shape. The equilibrium of this cylindrical shape is governed by the antago- 
nistic effects of elasticity and capillarity: elasticity tends to keep the film developable 
while capillarity tends to curve it in both directions so as to maximize the area of 
contact with the sphere. In the experiments, the contact area between the film and 
the sphere has cylindrical symmetry when the sphere radius is small, but destabilises 
to a non-symmetric, wrinkled configuration when the radius is larger than a critical 
value. We combine the Donnell equations for near-cylindrical shells to include a uni- 
lateral constraint with the impenetrable sphere, and the capillary forces acting along 
a moving edge. A non-linear solution describing the axisymmetric configuration of 
the film is derived. A linear stability analysis is then presented, which successfully 
captures the wrinkling instability, the symmetry of the unstable mode, the instability 
threshold and the critical wavelength. The motion of the free boundary at the edge 
of the region of contact, which has an effect on the instability, is treated without any 
approximation. 

1 Introduction 

The buckling of thin plates has been studied for a long time, both theoretically [U [2] 
and experimentally [3J. Initially, the main motivation was to avoid loads associated with 
catastrophic failure modes. Recent research efforts on thin plates and shells have been 
driven by technological applications involving thinner and thinner plates [4J, by the idea 
that controlled buckling can provide useful functionality [5], and by strongly non-linear 
phenomena appearing far above the bifurcation threshold [6j. The wrinkling of a semi- 
infinite elastic medium under finite compression, known as Biot's problem [8], as well 
as the buckling of a thin stiff film coated to a compliant substrate HU1 111) are just two 
examples of classical problems in mechanical engineering whose non-linear aspects have 
been well understood only recently. We refer the reader to |12j for a comprehensive review. 

When the adhesion between a thin film under residual compression and a thick sub- 
strate is relatively weak, buckling can take place along with delamination, resulting in the 
formation of blisters [131 H3J- Buckling-driven delamination can lead to various patterns 
which are affected both by the mode-mixity [15j . i. e. the dependence of the interfacial 
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Figure 1: Thin elastic films of thickness h, Young's modulus E and Poisson's ratio v are 
placed onto rigid spherical caps of radius p coated with a wetting liquid, with surface 
tension 7 = 22.4mN.m~ 1 . The liquid has been dyed in yellow to help visualization: the 
yellow region are the liquid meniscus. Note that the interior of the region of contact appear 
uncolored as there is almost no fluid there. Top view of three different experiments: (a) 
E = 2.8 GPa, h = 30 pm, p = 25 mm (b) E = 2.6 GPa, h = 15pm, p = 25 mm (c) 
E = 2.6 GPa, h = 15 pm, p = 60 mm. The parameter a, defined in equation (19), 
measures the strength of adhesion relative to the stiffness of the film, (d) Sketch of a cut 
through a vertical plane AA'. 



energy on the loading mode, and by the irreversibility of the interfacial fracture |16j . In 
recent experiments, a simpler variant of the classical delamination problem has been pro- 
posed, whereby the adhesion between the film and the substrate is provided by capillary 
forces arising from a thin liquid bridge \17\ [T8] : capillary forces are reversible and act like 
a self-healing interface crack. These experiments can be done at the centimeter scale. 

In the present paper, we study some buckling patterns produced by these experiments. 
Specifically, we consider the case of the capillary adhesion between a thin elastic film and 
a doubly-curved substrate. This geometry has been considered by one of us in a recent 
experimental paper [18J. It built up on previous work addressing the related case of a 
spherical shell adhering onto a planar substrate |19t [20l I21j . In our experiments, a rigid 
spherical cap is first coated by a wetting liquid and a thin polypropylene film is then 
applied onto it. As it wraps the sphere under the action of capillary forces, the film is 
forced to stretch by Gauss' theorema egregium |22j . Stretching allows it to make up for the 
mismatch of Gaussian curvature, which is zero in the planar film and non-zero along the 
spherical substrate. This leads to a variety of adhesion morphologies, as shown in figure [TJ 
where the contact region varies from a simple band to complex branched patterns. 

In reference [18] , one of us studied the antagonistic effects of elasticity and capillarity 
using order of magnitude arguments, and proposed an estimate for the size of the region 
of adhesion which successfully compares to the experiments. 

Here, we study these patterns quantitatively. In particular we address the transition 
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shown in the figure, whereby a band- like region of contact with straight edges (figure [T^i) 
bifurcates into a sinuous pattern with undulatory edges (figure [Tja) when the adhesion 
becomes stronger or the film becomes thinner. This is interpreted as a buckling bifurcation 
caused by compressive stress along the straight edges. We carry out a stability analysis 
based on the classical Donnell equations for nearly-cylindrical shells, modified to account 
for the effect of adhesion. The motion of the free boundary at the edge of the region of 
contact is considered without any approximation. 

The paper is organized as follows. In section [2j we derive the equations for a nearly- 
cylindrical elastic shell adhering to a sphere, with an emphasis on the equilibrium condi- 
tions along the edge of the moving contact region. In section [3j we derive a non-linear 
solution to these equations relevant to the unbuckled configuration with cylindrical sym- 
metry. These results are compared to experimental data. In section [4j we then study 
the linear stability of the cylindrical solution. The predictions regarding the symmetry of 
the buckling modes, their wavelength and the critical loads are compared to experimental 
data in section [H 

2 Governing equations: Donnell's shell equations with ad- 
hesion 

The Donnell-Mushtari-Vlasov equations for nearly-cylindrical shells, simply called the 
Donnell equations thereafter, are derived by combining the general equations for shells 
undergoing finite displacements with specific scaling assumptions for the magnitude of the 
displacement. Even though the Donnell equations have been known and used for a long 
time, this derivation is useful as it highlights the simple assumptions that underlie them. 
More importantly, the variational framework that is used to derive the Donnell equations 
allows us to include adhesion and the presence of a moving boundary in a natural and 
consistent manner. 

2.1 Geometry 

The reference configuration considered here is shown in figure[2^. A thin cylindrical shell 0, 
of half-length L x , half- width L y , thickness h and radius p rests tangentially to a sphere of 
radius p. The Lagrangian coordinates along the shell are denoted by x and y = p6, where 
8 is the azimuthal angle in the reference configuration. Let dfl^ and OQ^ correspond to 
the edges of the shell at x = ±L X and y = ±L y , respectively. We assume that the shell is 
made of an isotropic linear elastic material with Young's modulus E and Poisson's ratio 
v . 

A typical deformed configuration is shown in figure [2}). The displacement of the middle 
surface of the shell is denoted u x (x,y), u y (x,y) and w(x,y), where the first two functions 
are the tangential components of the displacement, and w(x, y) is the radial displacement 
(deflection), counted positive towards the exterior of the shell. In the deformed configu- 
ration, the position of a generic point lying on the center-surface of the shell reads 

x{x,y) = (x + u x (x,y))e x + (p + w(x,y))e 1 (9{x,y)), (la) 
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Figure 2: (a) The reference configuration: the cylindrical shell makes tangent contact with 
the sphere along the great circle x = 0. (b) A typical deformed configuration: the region 
of contact Ojy, shown in dark grey, is bounded by two curves G- and G + . 
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where the azimuthal angle in deformed configuration reads 
of \ y + u y (x,y) 

0(x,y) = y - . (lb) 



Here e 1 (^) = (0, cos 6, sin#) is the radial unit vector in the plane perpendicular to the axis 

The purely radial displacement w s < that brings a point from the cylindrical reference 
configuration onto the sphere is such that x 2 + (p + w s (x)) 2 = p 2 , and so 



W; 



:(X) = — (p — \f p 2 - X 2 ) 



where x is the axial coordinate and r is the radial coordinate. In the following, we shall 
use a second-order approximation, valid for \x\ <C p, 

w s {x) = -\* 2 -. (2) 
z p 

In this approximation, the sphere has been effectively replaced by its osculating paraboloid. 
The non-penetration condition is then expressed as a unilateral constraint, 

w(x) > w s (x). (3) 

We denote Qw the contact zone between the sphere and the shell, and Qp the free 
part of the shell, 0, = Qw U Neglecting the width of the meniscus, we assume that 
the boundary G = Q.w l~) between the two domains is made up of two curves, 

G = G-UG+. (4) 

In addition, we assume that each one of this curve is a graph, 

G-={(x,y)\x = g_(y)} (5) 
G+ = {(x,y) | x = g+{y)}. (6) 

This assumption is valid for configurations that are close to an axisymmetric configuration, 
for which both the functions g- and g + are constant. The wet region Q\y lies inside the 
boundaries G_ and G + and the free region Qp lies outside, 

n w = {{x,y) \ g-(x) <x < g + (y)} (7) 
= {(x,y) | x < g_(y) or g+(y) < x}. (8) 

In the rest of the paper, Greek indices such as a or f3 represent the Lagrangian coordinates, 
a = x or a = y, and follow the implicit summation convention for repeated indices. 

2.2 Scaling assumptions 

Before introducing the fundamental mechanical quantities, which are the membrane strain 
e a /3, the curvature strain b a p, the membrane stress n a( g and the bending moment m Qj g, we 
present the scaling assumptions that underlie the Donnell equations. 
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First, the thin-shell theory assumes that the following slenderness parameter is small, 



v 2 = t^- 0) 

yllp 

We shall therefore consider the limit 

77 <1. (10) 

The Donnell equations assume that the shell remains close to the cylindrical configura- 
tion of reference. How close depends on the small parameter 77: as shown in|Aj it is natural 
to rescale the deflection by prj 2 , in-plane lengths by prj, and in-plane displacements by 
pif . Therefore, we define rescaled lengths by 

x = — y = — , (11a) 

PV PV 

and the rescaled displacement by 

u a (x, y) = — ^ u a (x, y), (lib) 

w(x,y) = — Kw(x,y). (11c) 
pr] z 

The parabolic approximation of the sphere profile in equation ([2]) reads, in dimension- 
less variables, 

w s (x) = - l -x 2 . (12) 

Note that the assumption under which this approximation has been derived, |x| <C p is 
consistent with the new scales introduced here: when x is of order 1 (and in fact, as long 
as it remains smaller than the large number r/ _1 ), then |x| <C p and the sphere is indeed 
well approximated by its osculating paraboloid. 

For the sake of consistency, we rescale the membrane strain e a/ g by r? 2 , the curvature 
strain b a p by 1/p, the membrane stress n a p by C rj 2 and the bending moment m a /3 by 
B I p, where B = Eh 3 / [12(1 — u 2 )] and C = Eh/{l — v 2 ) denote the bending and stretching 
moduli of the shell, respectively. This is written 

e a p(x,y) = \ e a p(x,y) (13a) 

b a /s{x,y) = pb a/ 3(x,y) (13b) 
1 

m a p(x,y) = ^m a p{x,y). (13d) 

The definition of the small parameter rj given earlier in eq. ([9]) warrants that the bending 
energy density \ m a p b a fs and the stretching energy density | n a p e a p are commensurate 
when rj goes to zero. 



n aj 3{x, y) = n aP (x, y) (13c) 
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2.3 Membrane and curvature strains 



From now on, we shall use dimensionless quantities everywhere, and drop bars to easy 
legibility. The above scaling assumptions allow one to simplify the expressions for the 
membrane and curvatures strains from the general theory of shells as follows, 

e a p(x, y) = ( ^ La ^ x,y ^ > 1 u pA x >v) + ^ ^ w ^ + 1 w ^ Xf y ) w ^(x, y) (14a) 

b a p(x,y) = -w ta p(x,y). (14b) 

These expressions are shown to derive from the scaling assumptions in|Aj Here we use the 
Kronecker delta symbol, which is defined by 

8 aP = { 1 ifa = /? (15) 

and use commas in indices to denote partial derivatives, as in the expression 

du x (x,y) 



Ux, y (x,y) 



Oy 



The strain approximation in equation (14a) is at the heart of Donnell's model for 



nearly-cylindrical shells, and has been used by numerous authors, see for instance 
Only the non-linear terms that are important near the onset of buckling have been retained. 

2.4 Constitutive law and energy 

The shell's material is assumed to be linearly elastic and isotropic. The Hookean consti- 
tutive laws read: 

n a p(x,y) = (l-u) e a p(x,y) + v5 a/ 3 (trg(x,y)) (16a) 
m a p(x,y) = (1 - v)b a p{x,y) + v5 ali (trk(x,y)) (16b) 

This is a rescaled form of the original constitutive laws given in the Appendix in equa- 



tion (91). Thanks to our rescalings in equation (13), the stretching and bending moduli 
have effectively been set to one. 

In the experiments, the fim is naturally planar and is bent into a cylindrical shape 



by the action of capillary forces. By contrast, the constitutive equation (16b) describes 
a naturally cylindrical shell with natural radius p: the bending moment m a p is zero 
when the bending strain b cancels, which happens in the configuration of reference. The 
case of a naturally planar film could be treated by modifying bending strain, replacing 
b with (b + KQtiy ® e y ) everywhere, where kq = — 1/p is the curvature strain in reference 
configuration. As we shall see later, only the derivatives of the bending moment enter into 
the local equations of equilibrium. As a result, n$ is absent from these equations, and only 
affects boundary layers near the free edges. To simplify the analysis of these boundaries, 
we ignore this and set kq = 0. This is a valid approximation for the entire pattern, except 
very close to the film boundaries. 

The linear constitutive laws correspond to the following quadratic elastic energy, 

E shell = J J ^ (n>a/3 e a/ 3 + map b a p) dx dy (17) 
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which is the sum of a stretching and a bending term. 

The adhesion between the shell and the rigid sphere is modelled by an energy (27) 
per unit area of the contact region. Indeed the liquid perfectly wets both the shell and 
the rigid sphere, and any increase in the area of contact removes two air/liquid interfaces, 
each having an energy 7 per unit area |24j . where 7 is the surface tension of the liquid. 

The rescaled adhesion energy thus reads: 

E adhesion = ~ I I 2 a 2 xw(G;x, y) dx dy, (18) 

J JQ 

where Xw(G; x, y) is the characteristic function of the wet domain: xw(G; x, y) = 1 in the 
wet region (x,y) G and xw{G;x,y) = in the free region. The coefficient a 2 is the 
capillary energy rescaled by the typical energy per unit area introduced earlier, B/p 2 : 



B/li 1 \L 

Here we have introduced the elastocapillary length 




t -''--'' M [f-V (20) 

which arises from a balance of the bending rigidity of the shell B and the adhesion en- 
ergy 7. At length scales smaller than L ec , bending stiffness dominates capillary effects. 
The elastocapillary length sets the typical radius of the sphere beyond which adhesion is 
possible, as noted in reference [18] . 



2.5 Constraints and boundary conditions 

In the wet region, the condition of contact with the sphere reads: 

w(x,y)=w s (x) in Q w , (21) 
where w s is the parabolic approximation of the deflection of the rigid sphere derived in 



equation (12). We assume that the contact is frictionless. 

We consider that the shell is infinitely long in its y direction. This is captured by the 
following boundary condition on the remote edges dQ^ : 

u y (x,L±) = V ± (22) 

where the two unknown scalars V + and V~ denote an unknown rigid-body motion of 
the edges consistent with the symmetry. They will be determined later by a condition of 
equilibrium of the boundary. This boundary condition means that the remote edges dQ^ 
remain contained in a plane passing through the X axis, even though this plane can freely 
rotate about this axis to accommodate an average extension or contraction of the shell in 
its y direction. In the experiments, the film has a finite extent, and the above boundary 
condition is relevant to the case where its size is much larger than the wavelength of the 
instability: anticipating the notations of section [4j this writes L y 3> 2ir/k c . 

Finally, we introduce a new set of variables (q x , q y ) which is the local slope of the shell 
with respect to the mobile frame, 

Qa(x,y) = w a (x,y) (23) 
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In the following, these q a 's will be considered variables independent from the deflection w, 



and the equation (23) just written will be viewed as a constraint. This allows our second- 



order variational problem to be written as a first-order one having additional variables 
and constraints, and simplifies the derivation of the equilibrium equations. 

The Lagrangian of our constrained minimization problem is formed by augmenting 



the elastic and adhesion energies in equations (17) and (18) with the constraints in equa- 
tions (21), (22) and (23) by means of Lagrange multipliers, denoted n(x,y), v ± {x) and 



X a (x,y): 



E tot (u a ,w,q a ,V ± ,G,7r,is ± ,\ a ) = 

[[ 7; ( n a(5 e a p + m a/3 b a/ s) dx dy - [[ 2 a 2 xw(G) dx dy 

J JQ * J JQ 



// TT(w-w s )xw(G)dxdy-S~] 

J JQ ~i J d 



v s (u y - V s ) dx 



/ / K {w, a ~ q a )dxdy. 

J JQ 



(24) 



The interpretation of the Lagrange multipliers just introduced will be given in section [276] 
Note that the constraints (22) on the boundaries dVt~ and <90+ are taken care of using a 



summation over the values s = — and s 



2.6 Equations of equilibrium 

In this section, we derive the equilibrium equations and the boundary conditions using 
variational calculus, by cancelling the first variation of the Lagrangian just written. This 
yields the classical Donnell equations for shells inside each domain Vtw and Qp, as well as 
boundary conditions, including an adhesion condition at the interface G between the wet 
and free regions. 

We first compute the variation of the Lagrangian with respect to the unknowns (u a ,w, q a , V^, tt, u^, A c 
considering that domains £lw an d and their boundary G remain fixed, as expressed 
by the notation Sg = 0, 

8E\ ot = SE tot (5u a , 5w, 5q a , SV ± , 8g± = ^,5^,8^ ,8X a ) = 



II (n a p8e a p + 7n a pSb a p) dxdy- // ir xw (G) 5w dx dy 

JJQ J JQ 



E 



s=± J 9^y 



v s 5u y dx + 



v s dx 5V S 



(w - w s ) Xw{G) 8tt dx dy 



where we use the fact that 8{n a p e c 




X a {8w :a — 5q a ) dx dy 



V s ) 5 v s dx 



2n a p5e a/ 3 and 5(m al3 b al3 ) 



. w 



q a ) 8X a dx dy 



(25) 

2 m a p 8b a p, as we 



assume linear constitutive laws. The three last terms, which come from the variation 
with respect to the Lagrange multipliers (jr, v^, X a ), enforce the kinematic relations in 



equations (21), (22), (23), as expected 



The complementary variations with respect to the position of the boundary G will 
be computed later, with the help of [Bj They yield jump conditions across the boundary, 
which include the condition of adhesion. 



9 



Calculating the first variation of the strain defined in equation (14), and using the 
symmetry of the stress tensors n a p and m a p, we have n a p 5e a p = n a p (8u a! p + 5 ay 5p y 5w + 
wpdw^a) and m a p5b a /s = —m a p 8q a> p. Note that the symbol 5 without any subscript 
denotes the variation of a function, while 5 a p with subscripts is the Kronecker delta symbol 



introduced in equation (15). Integrating by parts and grouping the terms, we have: 



5E. 



lot 



J J ( [ n ap,p] 5u a + [it xw{G)+(n af3 w ja ) t p-n yy -\ aj a] 5w+ [-m a p t p-X a ] Sq^ 



dx 



+ / ( [n a i3-5 a y v s XdQ.* ] 8u a + [n a p w >a -Ap) 5w-[m a p] 5q a ) N/3d£+ 
Jan v £z± 7 
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where we denote the entire lateral boundary by Oft = <9f2^ U dQ^ , the element of length 
along the lateral boundary by d£ = dx on dQ^ and di = dy on dQ^ . In addition, 
N = (N x ,N y ) stands for the normal to a boundary <9Q, N being oriented towards the 
exterior of the domain 0,. 



The equilibrium condition, 5E. 



t 

tot 



0, yields the following equations on fi: 

n a /3,/3(x,y) = (27a) 

m a/ 3 >a /3(x, y) + n a p(x, y) w j(X p(x, y) - n yy (x, y) + xw(G; x, y) vr(x, y) = 0. (27b) 

Here we have replaced the shear force \ a by its expression X a = —m a p,p- The latter comes 
from the condition associated with perturbations 5q a . 

The equations (27) are the Donnell equations for thin cylindrical shells, see for in- 
stance [23]. The first equation (27a) is a tangential balance of forces. The second equa- 
tion (27b) is a transverse balance of force. The Lagrange multiplier ir can be interpreted 
as the pressure of contact with the sphere in the wet region. The equation (27b) allows 
one to compute the contact pressure ir(x, y) in the wet region where w = w s is known; in 
the free region, the pressure term ir(x, y) is zero and the equation is an equation for the 
unknown deflection w(x,y). 

On the lateral boundary dft£: we recover the natural boundary conditions for a plate 
or shell with a free edge, see for instance [25], 



m. 



n ax (±L x ,y) = 
m xx (±L x ,y) = 
m xy,y (±L x ,y)=0. 



(28a) 
(28b) 
(28c) 



Note that equation ( 28c ) comes from the integration by parts of 5q y in equation ( 26 ) , as 



= (Sw) tV on dfl^. 

On the other lateral boundaries d£ly, the boundary conditions read: 



n 



■.it/ 



(X,±Ly) 
l\ X 1 ^Ly) 



m 



n yy\ 

m yy (x, ±L y ) 

<<xy,x ^Ly) 



yy,y(.%i -^--^y) 2 

Thyy\X^ i-Z^y) dx 





v (x) 


— u (x)w t y(x, ±L y ) 
0. 



(29a) 
(29b) 
(29c) 
(29d) 

(29e) 
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Here we used again the fact 5q x = (Sw) jX on dQ,^. The Lagrange multiplier can be 



interpreted as the normal stress on dQy. The equation (29e) expresses the fact that no 



average force is applied on the shell in the y direction (natural boundary condition) 
2.7 Equations for the moving boundary 

We recall the general expression of the corner conditions in a two-dimension domain in 
[B] known as the Weierstrass-Erdmann conditions. They yield the jump conditions at 
the moving interface between two subdomains, in any minimization problem where the 
contributions to the objective function have different expressions in each subdomain. This 



applies to the Lagrangian of our problem in equation ( 24 ) , which can indeed be written as 



E to t= \\ C W {u a ,w,q a ,TT,X a )dxdy + // C F (u a , w, q a , n, A Q ) dx dy + / 
The integrands have different expressions in the free and wet regions, 



(30) 



C F (u a , w, q a , 7r, X a ; x, y) = - (n aj3 e a p + m a p b a p) - \ a (w }0l - q a ) (31) 
£ w (u a ,w,q a ,7r,\ Q ;x,y) = C F (u a ,w,q a ,ir, X Q ) - 2 a 2 - ir (w - w s ). (32) 



These notations conform with those of equation (96) in [Bj As usual with shell models 
having non-zero bending rigidity, the tangent displacement u a (x,y) is required to be con- 
tinuous, and the transverse displacement w(x,y) to be C 1 -smooth. In particular, across 
the boundary G, we have 

M = (33a) 
M = (33b) 
HJ = 0, (33c) 

where 

lA} = A F (xG,y G )-A w (x G ,y G ) (34) 

denotes the discontinuity of a function A across a point (xcVg) ly m g on the boundary 
G. 



By differentiation of the equalities (33) along the boundary G, we find 



KtI = (35a) 

Karl = 0, (35b) 

where T = (T x ,T y ) is the unit tangent to G, and a T in subscript denotes the tangent 
derivative fx = T a f a . The functions w and 

Qa — have prescribed values in the wet 
region £lw but can take any value in the free region tip. By contrast, u a is unconstrained 
in the entire domain 0. 



Using the fact that the elastic energy is a quadratic form of the strain by equation ( 17 ) 



first, and using the definition of the strain in equation (14) next, one can compute the 



so-called generalized momentum and identify the result with the internal stress, up to a 
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sign, 



3D 



de 



n 



IP 



1 

■ - fP j 



db 



dq 



m. 



ip 



a,/3 



IP 



dq, 



Sap Spy) 



-m 



a/3 



(36) 
(37) 



where % = W, F is any of the wet (W) or free (F) region. 

We now apply the corner conditions derived in[B]after identifying the regions Oi = £lw, 
0,2 = &f- The unknowns, collectively denoted £ a in the appendix, are the in-plane 
displacement up, the deflection w and the slope qp. 



When the equilibrium condition (101) derived in the appendix is applied to the un- 
known ia = Ua-, we find \n a p\ Np = \n a N\ = after using equation (36). Note that Np 
refers to the local normal vector, while n a p refers to a generic component of the membrane 
stress. For a = T and a = N this yields [n^] = an d [^Afiv] = 0, which are the classical 
conditions for the in-plane equilibrium of the boundary. 

It turns out that the third independent component utt of the membrane stress is also 
continuous across the boundary, even though this does not directly follow from equilibrium. 
To show this, let us write the constitutive law (16a) in dimensionless form, utt = vriNN + 
(1 — v 2 ) 

ctt- In the right-hand side, unn is continuous as we have just shown, while &tt 



is continuous as a consequence of the smoothness conditions (33) and (35). Therefore, the 



tangential stress utt is continuous as well. To sum up, we have shown that all components 
of the membrane stress are continuous, 



\n a p\ = 0. 



(38) 



Using the constitutive law for stretching in its inverted form, one shows that the membrane 
strain is continuous as well, 

Ml = 0. (39) 

As a result, the density of stretching energy is continuous, n a p e a p] = 0. 

The discontinuity in the energy density is therefore the sum of an adhesion term, which 
is present only in the wet part, and a discontinuity in bending energy, 



[£]= C w = 2a 2 + ^{m a pb, 



1 



(40) 



Considering now the condition (102) for equilibrium of the boundary with respect to 
the variable £ a = qa = w,a> we have 



1 



2 a 2 + - \m a p b a pj - m*p {b^j Np Np, = 0. 



(41) 



It is consistent to use this equation (|102[) as the value of q a = w y a = w s a is prescribed in 
the adhering region. 

Upon insertion of the constitutive law first, and of the continuity relations for the 



curvature next, one can rewrite the second term of equation (41) as 
\ \m a p M = \ {b 2 NN + b\ T + 2 v b TT b NN + 2 (1 - u) b 2 TN \ = X - ({b 2 NN j + 2vb TT [6 



NNj) ■ 

(42) 
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Note that only the AW-component of b a p = —w >a p is possibly discontinuous. As a result 
only the set of indices 7 = N and f3' = N need be considered in the last term of (41) 



Using the constitutive law for bending ( 16b ) again, we can rewrite this last term as follows, 
- m F p 1 M Np N p , = -{b F NN + v b TT ) [b NN j . (43) 



Inserting the expressions \A2\ and (43) into the equilibrium condition (41), we find 



that the terms proportional to Poisson's ratio cancel out: 



1 



2^ 2 + {{b F NN ?-K N ?)-b F NN (b 



A N ~ ^NN 



)=0. 



(44) 



After using the definition of the jump operator in equation (34), this can be simplified as 

1 



2<r*-^[b NN ] 



0. 



(45) 



By developing in the non-penetration condition in equation ([3]) in Taylor series up to 
second order, and by using the continuity of the deflection and the slope in equations (33), 
we have 

w,nn(xg, Vg) > w StNN (x G , y G ). (46) 

In terms of curvature, this yields b F JN (xG, Vg) < b^ N (xG, Vg), which we rewrite using the 
jump notation as 

{bNNj < 0. (47) 



Thus, the equation (45) implies 



lb 



NNtt 



-2 a, 



(48) 



as a is a positive number. 

To sum up, there are seven independent continuity and jump relations that must be 
enforced at the moving interface g. They read 



["^ Ja;=g(j/) 
\Uy\ x =g(y) 
Mx=g(y) 

lw,N}x=g(y) 

\nTN\x=g{y) 
l n NN\x=g{y) 
lmNNjx=g(y) 












-2a. 



(49a) 
(49b) 
(49c) 
(49d) 
(49e) 
(49f) 
(49g) 



Equation (48) comes from equation (49g) after using the constitutive law. The notation 
x = g{y) emphasises the fact that the jump operator is defined on the boundary G whose 
equation is x = g±(y). A term proportional to the perturbation 5g± of the boundary will 
appear in these equations when we consider the linear stability later on. This perturbation 
8g± must be considered as the shape of the boundary may be affected by the buckling. It 
corresponds to the quantity a\ introduced later. 



The jump relation (49g) arising from adhesion has previously been derived in an ef- 

, [27] , as well as for elastic plates [131 [2S] ■ This equa- 



fectively ID (axisymmetric) case 



tion (49g) is also known in the context of classical fracture mechanics, as the delamination 
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of a thin film on a rigid substrate can be seen as the propagation of a crack along an 
interface [15]: the energy release rate G per unit width of a film adhering on a curved 

substrate reads G = ■ 
G = 2 7 . 



and equation (49g) expresses from the balance of energy, 



3 A non-linear solution for the unbuckled state 

In this section we derive an axisymmetric solution to the equations derived in section [2j 
The wet region Ovp is then a strip bounded by two symmetric circles with equation x = 
9±{v)i where g±(y) = ±cto- Here ao denotes half the width of the strip, and will determined 
later as a function of the adhesion number a. 

In the axisymmetric case, the displacement is of the form 

u x (x,y) = u° x (x) (50a) 
u y {x,y) = E y (50b) 
w(x,y) = w % (x), (50c) 

where the index i denotes the region (i = W if |sc| < ao and i = F if \x\ > ao), the symbol 
'0' refers to the unbuckled axisymmetric state. The unknown constant Eq measures the 
uniform hoop strain, and will be determined later. It is related to the unknown uniform 
tangent displacement at the edges y = ±L y by Eq = u y>y = (V + — V~) / '(2 L y ). 
Imposing a linear dependence of u y on the azimuthal variable y warrants that the hoop 
strain will be independent of y, as required by the symmetry. 
We compute the membrane strains 

e° xx (x, y) = uf(x) + - K'(x)) 2 (51a) 

e° xy (x,y) = (51b) 
e° v (x,y)=E + v%(x). (51c) 

They are all independent of the coordinate y, as required by the symmetry. All bending 
strains are zero, except for the axial component 

k°(x,y) = b xx e x ®e x , b° xx = -w l "(x). (51d) 

By the constitutive law for stretching we have n xy (x,y) = and n a p jy (x,y) = 0, as 



required by the symmetry. The equilibrium condition (27a) along the direction x (a = x 



in this equation) implies that the stress n xx does not depend on x either. By the boundary 



condition (28a), this quantity vanishes everywhere, n xx (x,y) = 0. In terms of the in-plane 
strain, this writes e xx + v e yy = 0, an equation which will be useful to reconstruct the 
in-plane displacement u x l (x): 

uf(x) = -^K'(x)) 2 - vwi(x) - uE . (52) 



2 



Elimination of e xx = —ve yy from the constitutive relation (16a) yields the expression of 
the only non-zero stress component, 

U {x,y)=n yy e y ®e y , n° yy = (1 - v 2 ) e° yy (x, y) = (1 - v 2 ) (E + w l (x)). (53) 
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In the wet region \x\ < an, the contact condition reads 



1 



(54) 



The membrane stress there is then found by inserting this expression into equation ( |53[ ), 
up to the constant Eq that will be determined later. The in-plane displacement u° (x) 
is found by integration of equation (52) with the initial condition u x (0) = imposed by 
the symmetry. 

By symmetry, many term cancel in the equation (27b) for the transverse equilibrium 
in the free region \x\ > qq, which reads 



1 



1 



F"" t \ i Ft \ 
w (x) + w [X) 



-En. 



(55) 



We consider the case of a shell of infinite length in the x-direction, L x — > oo. This 
corresponds to the situation in the experiments where the length L x of the film is much 
larger than the width of the wet region, that is much larger than the typical in-plane length 
\fph which we used to make lengths dimensionless. The generic solution of equation (55) 
that is bounded near x — > ±oo reads 



-Eq + e ( A\ cos 



where A\ and A2 are unknown amplitudes, and 

4 



x*{v) 



+ A2 sin 



x*(u) 



1 



1/4 



(56) 



(57) 



is a scaling factor applicable to in-plane lengths. This x* (v) is a known function of Poisson's 
ratio. 

In order to determine the four constants (A\, A2, Eq, qq), we use all the boundary 



conditions that are not automatically satisfied, namely the continuity conditions (49c 



and ( 49d ) , the equilibrium condition for the average force along the y-direction ( 29e ) , and 
the jump condition (49g) depending on the adhesion number a: 

w o( a o) = w^iao) 
Oo) 

n [ ' (x) dx = 



F' 1 \ W'/ 
w (oo) = w 



(58a) 
(58b) 

(58c) 

(58d) 



Even though some of the equations for the problem are non-linear, equations (58) happen 
to be linear with respect to the unknowns A\, A2 and Eq, when we make use of equa- 
tions (53), (54) and (56). This leads to a set of linear equations whose coefficients depend 

(59a) 



non-linear ly on their arguments, 



M 



(1Q 



,<t ) -X = 



where 



X= (A l ,A2,Eo,ao i 



(59b) 
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and 



A/[S= * a 

x*{u) 



I e~ a cosd e~ a sin a — 1 s \ 

-e~ a (cosd + sind) e _a (cosd — sind) 1 

e~ a (cosd — sind) e _a (cosd + sind) 

\ e _a sind — e~ a cosd ^ (| + cr) / 



a 
_ 'A 
3 



A necessary condition for this linear system to have a solution is 



(59c) 



n=UHh a (60) 

One can simplify the determinant of M and rewrite this equation as 
3 + 6d + 6d 2 + 2d 3 , „ a 

= <r, where a = — - (61 

6(1 + a) x*(y) 

This implicit relation, plotted in fig. [3^, selects the half-width ao = (x*(y)a) of the 
region of contact Q w in the axisymmetric configuration, as a function of the dimensionless 
adhesion number a. Note that this half-width ao goes to zero as a decreases to the value 
1/2: the adhering solution disappears when the adhesion is too weak, a < 1/2. 

From now on, we assume a > 1/2 and consider the value of ao that is the unique 



solution to equation (61 ). Then, the matrix in equation (59a) is singular and the solutions 



X of the linear system span a line. Generically, there is a unique vector X_ whose last 



component equals the square of the quantity ao just found, as imposed by equation (59b). 
The other components of this particular vector X_ set the values of the unknowns A\, 
Ai and Eq. For any value of the adhesion parameter a > 1/2 this defines a unique 
axisymmetric solution. Some of these solutions are represented in figure^, for particular 
values of a. The residual stress n°(x) given by equation (53) is plotted in figures ^> and 



c. We note that this residual hoop stress is compressive near the edge of the region of 
contact, as the film is pulled towards the axis by the adhesion. 

In figure HI we compare the prediction for the width of the contact region in equa- 



tion (61) to the experimental data taken from Ref. [18], with no adjustable parameter. 
We find a good agreement. Interestingly, the agreement is good even in the post-buckled 
regime (when a is larger than approximately 5), if we measure ao in the experiments as 
half the average strip width; this is unexpected as the present solution is not applicable 
above the bifurcation threshold. 

In the limit of strong adhesion, when a is large, the half-width ao predicted by equa- 



tion (61) becomes large as well, and is given asymptotically by ao = v3cra;*(^). Restoring 



the physical units, this yields 

■^Mikr- < 62 > 

We recover the scaling form proposed and verified in reference [18] , and have obtained the 
value of the coefficient in addition. 

We have derived a family of non-linear solution of the Donnell equations analytically 
that describe unbuckled, axisymmetric shapes of the film. The equations are non-linear 
and it is remarkable that these solutions can be derived without approximation. There is 
a unique axisymmetric solution when the adhesion is large enough, a > 1/2. For lower 
adhesion, a < 1/2 no axisymmetric solution with extended contact can be found. 
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Figure 3: Axisymmetric solution, (a) Half-width of the band-like region of contact, as a 
function of the reduced adhesion a. The vertical axis on the left-hand side shows the value 
of aQ/x*(v) which is independent of Poisson's ratio, and that on the right-hand side shows 
the corresponding value of ao for the particular case v = 0.4, as in the experiments, (b) 
Hoop stress n° (x)/n*(v), rescaled using n*(u) = l/(x*(u)) 2 , as a function of the rescaled 
axial coordinate x/x*{u), for increasing values of the adhesion a. When expressed in these 
variables, the curves are independent of v. Dashed lines correspond to the wet region, solid 
lines to the free region. The different curves correspond to different values of adhesion, 
namely a = 1,3,5,7 and 9. (c) Profile of the shell as obtained by cutting through a 
plane passing through the axis, for increasing values of the adhesion number a. Only one 
half of the cut is shown, the entire profile being symmetric with respect to x = 0. The 
profiles corresponding to different values of a have been offset vertically for clarity. As 
the adhesion a increases, the width of the region of contact increases, and the hoop stress 
becomes more and more compressive around the edge of the region of contact. This points 
to the existence of an instability, which is studied next. 
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100 



Figure 4: Size of the contact zone as a function of the adhesion parameter a. Circular 
symbols correspond to experimental data from Ref. [18]. The solid line corresponds to 



the prediction of equation (61), and the dashed line to the asymptotic behaviour of equa- 



tion (61). As the solid curve comes from the analysis on the unbuckled stated, it cannot 
be expected to be accurate far above threshold, when a becomes larger than ~ 10; a good 
agreement is observed nevertheless. 



4 Stability analysis 

Given the presence of compressive hoop stress in the axisymmetric solution shown in 
figure |3]3, the neighborhood of the edges of the region of contact can become unstable, 
especially when the adhesion number a becomes large. Making use of the explicit solution 
for the unbuckled state, one can approach this question by studying the linear stability of 
the unbuckled state. This is the goal of the present section. 



4.1 Perturbations 

We investigate the presence of bifurcated branches near the unbuckled state, and introduce 
a perturbation of the previous solution in the form: 

u x (x,y) = ■u° i (x)+n^(x) cos(ky) (63a) 
u y (x,y) = E y +u 1 y l (x) sin (k y) (63b) 
w(x,y) = w l (x) +w\{x) cos (ky). (63c) 

Here the index 1 refers to perturbations from the axisymmetric state. As the axisymmetric 
solution is invariant in the y direction, the harmonic dependence of the perturbations 
on the y variable introduced above is the only one that we need to consider: a generic 
perturbation can be recovered by linear superposition. 

The edge G of the wet region may deform upon the instability. Therefore we introduce 
a perturbation of its boundary, 

9+{y) = ao + a x cos(ky). (64) 

Since the base solution is mirror-symmetric with respect to the plane x = 0, we shall only 
need to consider perturbations that are either symmetric, or antisymmetric. The benefit 
is that we only need to solve the linearized equations on half the domain, x > 0, using 
initial conditions at the center of symmetry x = that reflect the type of symmetry under 
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consideration — details will be provided below. The shape of the boundary at x = —ao can 



be reconstructed from the other boundary using equation (64): g~(y) = —ao — a\ cos(ky) 
in the symmetric case (also called the varicose pattern, as the width of the wet region gets 
modulated while its center- line remains straight), and by g+(y) = —ao + a\ cos(ky) in the 
antisymmetric case (in this sinuous mode, the width of the wet region remains constant 
to first order but its center- line undulates laterally), see figure [5] 

4.2 Linearized equilibrium in the interiors of the domains 



Inserting the expansions ( |63a ]63c ) into the equations of equilibrium (27), and linearizing 



around the unbuckled state (u x l , Eq, Wq, ao), yields three coupled, linear ordinary differ- 
ential equations for the functions u^(x), Uy(x), w\(x). These equations are fourth order 
with respect to w\ and second order with respect to and u y l . 

In the wet region the deflection is prescribed by equation ^ and (x) = 0. 
We are not interested in computing the perturbation to the contact pressure, and will not 



use the transverse equilibrium (27b) there: we are left with the linearized equations for 
in-plane equilibrium, which are two second order, ordinary differential equations for u x l 
and u y l . 

The linearized equations of equilibrium can be cast into an equivalent first-order form 
by introducing the state vectors, defined in each region by 

U w (x) = [ul w ,ul w ',ul W ,ul w '] (65a) 

tt ( \ r IF IF' IF IF' F F' F" F"'i 

y_p[x) = [u x ,u x ,u y ,u y , w{ , w 1 , w 1 ,w 1 \. (oobj 
The linearized equilibrium in the interior of each domain then reads 

U!w( x ) = Aw^ v ' a > x ) ' —w ( x ) for x e ^ w ( 66a ) 

V_f{ x ) = Ap( v i a > k, x) • U_f{ x ) f° r x ^ ^f, (66b) 

where A w (respectively A p ) is a 4 x 4 (respectively 8x8) matrix whose coefficients depend 
on the dimensionless coordinate x, on the adhesion parameter a, on Poisson's ratio v, as 
well as on the wave number k. This dependence arises either directly, or indirectly through 
the unbuckled solution (u° 4 , Eq, Wq, ao). 

4.3 Asymptotic behavior far from the sphere 

We first integrate the linearized equations of equilibrium along the free region, ao < x < 
L x . We start from the free end x = L x where we use initial conditions consistent with the 
stressfree boundary conditions, and proceed towards the moving interface x = ao- The 
perturbation computed near x = ao will be ultimately reconciled with that coming from 
the wet region using the equations at the mobile interface. 

We consider an infinitely long shell, L x — > +oo. This is an accurate approximation 
as in the experiments the film is much wider than the region of contact. We must only 
consider solutions of the linearized problem that remain bounded for large x: we start 
by studying their asymptotic behavior, which generically is exponential. The trial form 
u x F (x) = U X F exp (tx), u y F (x) = Uy F exp (rx) and wf(x) = Wf exp (rx) is inserted into 
the linearized equilibrium (66b). Denoting = [U X F , Uy F , Wf] the vector collecting the 
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unknown amplitudes, the equations of equilibrium are asymptotically satisfied provided 
the following condition holds: 



AJv,a,k,T)-U^=0, 
where A captures the asymptotic form of the linearized equations, 



(67) 



'^(-1 + v)+t 2 |(l + ^)r 

-f(l + i/)r 

— VT 



-k 2 



H-1 + v)t 2 



-k 



■1 



VT 

-k 

k A -l-k A + 2k 2 T 2 - r 4 
(68) 



The acceptable values of the decay rate r are found by requiring that equation (67) has 
non-trivial solutions: 



detA =-(-l + v)(k i 
=oo 2 



4fcV + r 4 + 6Fr 4 - v 2 t a - 4k 2 T b + r 8 ) = 0. (69) 



Each one of the eight complex roots, denoted Tj, is associated with an eigenvector U_ 



Fi 



The stressfree boundary conditions (28) will be automatically satisfied, provided the 



perturbation stays bounded for large x. Therefore, out of the eight possible exponential 
behaviors, we can keep all four non-divergent solutions which are such that the real part 
of Tj is negative, 3?(tj) < 0. By convention, these values of Tj are indexed by 1 < j < 4. A 



generic, non-divergent solution of the linearized equilibrium (66b) is then found by linear 
superposition, 



(70) 



an approximation that is accurate for large values of x m . The four unknown amplitudes 
of the converging modes are collected into a vector I F = [ii, 12,13, 14]. 



We can now use the asymptotic form (70) as an initial value, and integrate the lin- 



earized equation towards the edge of the free region, x = do- By linearity, the state vector 
C/p(a^~) depends linearly on the asymptotic amplitudes I F : 



HF( a o ) = §. F ( U > a > k ) " If- 



(71) 



This 5„ is the so-called shoot matrix S„, and its size is 8 x 4. Its columns are computed 
as follows. First the initial condition for the linearized equilibrium (66b) are computed 



based on the asymptotic form of the perturbation ( 70 ) , using a finite but numerically large 



value of x m ; all the coefficients ij's are set to zero except for one which is set to one. The 
linearized equilibrium is then integrated from x = x m to x = ao numerically, and the final 
state vector U F (ao) is used to fill in the corresponding column of S_ p . This shoot matrix 
captures the linearized response of the free region, including the remote stressfree edge, 
as viewed from the edge G. 



4.4 Symmetry conditions at the center of the domain 

As explained earlier, the symmetry of the base solution by a mirror reflection changing x 
to (—x) allows us to consider perturbations that are either symmetric or antisymmetric, 
while retaining full generality. These two types of buckling modes are depicted in figure [5j 
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(a) (b) 

Figure 5: Two types of buckling modes: (a) Symmetric (varicose) and (b) antisymmetric 
(sinuous) perturbations. 

The type of symmetry dictates the initial condition at the center of symmetry, x = 0. A 
varicose perturbation, shown in figure^, is denoted by a superscript (+) as it is symmetric 
with respect to the x axis; for this type of symmetry, the conditions u^ w (0) = Uy W (0) = 
hold. A sinuous perturbation, shown in figure [SJ), is denoted by a superscript (— ) as it is 
antisymmetric with respect to the x axis; for this other type of symmetry, the conditions 
u x W (0) = u y (0) = hold. We can thus define two independent initial state vectors 
in the wet region as follows: for the analysis of the varicose mode, f/^ 1 (0) = [0, 1,0,0] 
and U^} 2 (0) = [0,0,1,0]. For the analysis of the sinuous mode, t/^ 1 (0) = [1,0,0,0] and 
U^} 2 (0) = [0, 0, 0,1]. A generic initial condition compatible with the symmetry is obtained 
by linear superposition, using two unknown amplitudes which we denote is and i%: 

2 

^(o) = &+^'(o). ( 72 ) 



Integrating equation (66a) across the wet region Qp using each of the various modes in 



equation (72) successively, we compute a shoot matrix for the wet region: 



U^)=§^(p,a,k)-I w . (73) 

There are in fact two such shoot matrices, one for each type of symmetry, and with size 
4x2 each. 

4.5 Assembly 

Let us define an assembled state vector at the boundary x = ao by collecting those relevant 
to the free and wet regions: 

U(ao) = [U F (4),uU%),ai}- (74) 

We have also appended the perturbation a\ of the boundary shape, which will soon be 
needed. 

We can similarly define an assembled shoot vector by 

L=\L F ,L W M- (75) 
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This / is the main unknown of our stability problem, and will be shown to satisfy an 
eigenvalue problem. 



Then equations (71) and (73) can be rewritten in compact form as: 



(76) 



where the global shoot matrix Sj^ F is formed by assembling into blocks the shoot matrices 
previously computed: 



Sv 
1/ 



(77) 



Its size is 13 x 7. 



4.6 Equilibrium of the moving boundary 

At this point, we are ready to close the formulation of the stability problem by using 
the remaining continuity and jump conditions at the edge of the contact region. Upon 



linearization, the seven conditions (49) can be written in matrix notation as: 



g(v, a, k) ■ U{a ) = 0. (78) 
The matrix C is of size 7x7. It collects the coefficients appearing in these linearized 



equations. The equations (49) hold on the mobile curve x = g+(x) which is perturbed 



according to equation (64). As a result, the linearized equations have terms proportional 



to the boundary perturbation a\ times the gradients of the base solution, and these terms 
are used to fill the last column of C. This allows the perturbation to the boundary to be 
treated without approximation. 

Note that the inhomogeneous term (—2 a) appearing in the right-hand side of the 



adhesion condition (49g) disappears upon linearization. 



4.7 Linear stability formulated as an eigenproblem 



By combining the equilibrium of the interface in equation ( 78 ) and the integration in the 



free and wet domain captured in equation (76), we can write the linear stability problem 
as an eigenvalue problem, 



\g(u,a,k)-S^ F (u,a,k)] 1 = 0. 
The existence of a linearly unstable mode requires 

det \C(u, a, k) • §_^ F {v, k)] = 0. 



(79) 



(80) 



This is an implicit equation for the wave number k as a function of the adhesion parameter 
a and of Poisson's ratio v. The corresponding eigenvector I allows the linearly unstable 
mode to be reconstructed. 
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Figure 6: (a) Linear stability diagram for v = 0.4. The most unstable symmetric and 
antisymmetric modes are defined by a^r and k^r (b,c) Dependence of the critical wavevector 



kf and adhesion number a^r on Poisson's ratio. 



4.8 Numerical results 

We explained how the shoot matrix S^^u, a, k) and the linearized equilibrium of the 
boundary C(i/, a, k) can be computed numerically for specific values of their arguments. 
For a given value of Poisson's ratio v and for each type of symmetry, we repeatedly 
computed the determinant for different values of a and k, and then plotted the implicit 
curve defined by equation (80) in the plane (k,a). The result is shown in figure [6^, for 
a typical value of Poisson's ratio, v = 0.4. Numerical convergence with respect to the 
parameter x m required in the calculation of the shoot matrix S p was attested by the fact 
that further increase of x m did not significantly affect the stability curves; typically, this 
required x m ~ 12. 

The most unstable wavenumber k^ and the critical adhesion correspond to the 
minimum of each of the curves in the plane (k, a). For each type of symmetry, the most 
unstable mode depends only on Poisson's ratio, hence the notations c^(z^) and k^r(v). 
The corresponding curves are plotted in figure [6jo and c. 

The stability analysis predicts that the most unstable mode is always the antisymmetric 
(sinuous) one. Indeed, the curve corresponding to this mode is always below the other 
curve in figure [6J3. This is consistent with the experiments, where the varicose mode has 
never been observed. 

It turns out that the dependence of both a c ±(v) and k^{v) on Poisson's ratio can be 
captured by simple formulas which match the numerical results perfectly, within numerical 
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accuracy: 



0.718 
0.891 



<r~(y) ps 3.788, 
cr+(^) 4.184, 



(81a) 
(81b) 



where x*(z^) is the function defined in equation (57). We have no explanation to offer 
for the fact that cr^f is independent of v, and that k^(v) depends on v as l/x*(y). This 
probably points to the fact that a proper rescaling of the various quantities would allow 
the parameter v to be removed from the linearized equations altogether. We have not 
been able, however, to identify such a rescaling. 

The predictions relevant to the experiments, where we used a film with Poisson's ratio 
v = 0.4, are 



a' (OA) = 3.788, K(0A) = 0.486, 



2vr 



k c - (0.4) 



12.93. 



(82) 



Here i~ denotes the wavelength of the sinuous instability. These results are in qualitative 
agreement with the picture shown in figure [TJ the straight contact region looses stability 
somewhere between a = 1.37 and a = 4.01, and the buckled pattern is indeed sinuous 
(antisymmetric). A detailed and systematic comparison to the experiments is presented 
in the following section, using a novel experimental setup that allows the adhesion number 
a to be varied continuously. 



5 Comparison to experiments 

In order to check the predictions of the linear stability analysis, we have developed a novel 
experimental setup that allows the adhesion parameter a to be continuously varied while 
the shape of the contact area is monitored. To do so, we replaced the rigid spheres used in 
previous work [TBJ by an inflatable latex membrane. The membrane was cut out in a latex 
sheet with Young's modulus E m = 1.2 MPa and thickness h m = 0.6 mm, and glued on top 
of vertical rigid cylinders with radius p cy \ = 20, 38, 100, 295 mm, see figure [7k. Increasing 
the pressure by Ap inside the cylindrical container allows the radius of the latex sheet to 
be varied continuously from p = oo when the membrane is flat, to p = p cy \ when it is half 



a sphere; the adhesion parameter a then varies according to equation (19). We checked 
that the shape of the membrane is almost spherical. The value of the radius of curvature 
p(Ap) was measured from pictures taken from the side. 

The spherical cap was then coated with ethanol, whose surface tension is 7 = 22.4 mN m" 
Thin polypropylene films were applied on top of it, as shown in figure [Tjo. The films, pro- 
duced by Innovia films, have a Young's modulus E = 2.6 ± 0.2 GPa, a Poisson's ratio 
v = 0.4, and their thickness ranges from h = 15 pm to 90 pm. Ethanol was chosen because 
it wets both latex and polypropylene, and because its surface tension is not very sensitive 
to impurities. Snapshots are simultaneously taken from top and from the side to monitor 
the shape of the pattern as a function of the radius of the spherical cap. The transition 
from band-like contact pattern to sinuous pattern was observed as a result of membrane 
deflation, see figure [7[3,d. 

The typical evolution of the edges profile as a function of the radius p is shown in 
figure [8k. The edges remain straight until the sphere reaches a critical radius p c and then 
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(a) 




(b) 






Figure 7: (a,b) Experimental setup. A square latex sheet glued on top of a rigid cylinder is 
inflated by a pressure Ap and deforms into a spherical cap. It is coated with ethanol and a 
thin film is applied onto it. (c,d) Experimental images of the adhering regions, as viewed 
from top: the darker strip corresponds to the contact zone. Decreasing the pressure Ap 
increases the adhesion parameter a, and induces a transition from straight to oscillatory 
edges. 
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Figure 8: (a) Evolution of the edges as a function of the radius of the spherical cap p. All 
curves correspond to the same film, h = 50 /im. (b,c) Evolution of the rescaled amplitude 
ai of the edge profile as a function of the reduced adhesion number a. Variations of a 
have been achieved by deflating the membrane. The experiment has been repeated three 
times using the same film with thickness h = 50 /im, and the symbols correspond to the 
different experiments. The value of the threshold is very reproducible, but the amplitude 
of buckling in the post-buckled regime is more scattered. 
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Figure 9: (a) Critical radius of sphere p c leading to undulations of the contact region, as 
a function of the elastocapillary length L ec . Each data point corresponds to a specitif film 
thickness, and is obtained by detecting the threshold of instability in the experimental 
curve for a\(p), as in figure 03. The various data points correspond to thicknesses h 
ranging from 15 /im to 90 fim. Note the collapse of these points onto a single curve, (b) 
Wavelength of the oscillations £ c = 2i:/k c at threshold (p = p c ). In (a) and (b), the solid 
lines show to the prediction of the analysis of linear stability in equation ( |83| ) for v = 0.4, 
with no adjustable parameter. 



become undulatory; the amplitude of undulation increases as p is further increased. The 
observed buckled patterns are always sinuous, as predicted by the stability analysis. By 
fitting the edges profile with a cosine function a\ cos (kx + 4>), we extracted the wavelength 
and the amplitude of the oscillations. Our measurements for a particular polypropylene 
film (h = 50 pm) are collected in figure J8j}, which shows the dependence of the amplitude 
of oscillation on the sphere radius. The amplitude of the undulations goes to zero at 
the instability threshold, which is consistent with a supercritical bifurcation. When the 
experiment is repeated several times using the same film, the value of the bifurcation 
threshold is found to be quite reproducible. The scattering of the data for the buckling 
amplitude well above the threshold may attributed to variations in the volume of ethanol 
used, to friction between the film and the sphere, and to variations in the way that the 
film is initially laid out on the sphere. 

Restoring the physical units, we rewrite the predictions of the analysis of linear stability 



in equation ( 82 ) as 

p c (0.4) = 3.788 L ec , 4(0.4) = 13.52 a/ h L ec . (83) 

We recall that p c and £ c denote the radius of curvature of the sphere, and the wavelength 
of the pattern at threshold, respectively. 

In figure[9]we compare these predictions to the experiments. The elastocapillary length 



L ec is computed from the experimental parameters using equation (20). The wavelength 
of the pattern agrees very well with the prediction of the analysis of linear stability; there 
is no adjustable parameter. The agreement concerning the instability threshold is not as 
good but satisfactory; as apparent in Figure [7|1, the actual threshold is lower by about 
15% than that predicted by theory. This can probably be attributed to imperfections in 
the planarity of the film and in the way that it is layout onto the substrate, as well as to 
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the finite size of the meniscus. 



6 Conclusion 

We have investigated the adhesion of a thin film on a spherical substrate. The Donnell 
theory of nearly cylindrical shells has been used, taking into account the energy of adhesion 
with the sphere. We have derived the boundary conditions holding at the edge of the mov- 
ing boundary between the wet and free regions. A family of nonlinear solutions describing 
the unbuckled configuration of the film has been constructed. These solutions, which have 
cylindrical symmetry, are indexed by a dimensionless adhesion number a = p/L ec , that 
compares the radius of curvature of the sphere p to the elastocapillary length L ec . As the 
adhesion number a increases, a compressive stress builds up along the edge of the contact 
zone and makes the film buckle. We have carried out a linear stability analysis of the 
unbuckled solution, taking into account the motion of the moving boundary between the 
wet and free regions, and found of a critical value of the adhesion number a above which 
the contact region becomes sinous. The symmetry of the unstable mode, as well as the 
instability threshold and the wavelength found by the theory are in good agreement with 
the experiments. 

In our system, buckling is driven by the geometric frustration due to the mismatch 
between the Gauss curvature of the substrate and that of the planar film; a related type 
of buckling driven by geometric frustration has been investigated in plates undergoing 
swelling [29J, and in annular origami models having a curved crease [30 . Our paper 
provides a detailed analysis of the initial buckling of our system. Far above the buckling 
threshold, the buckling patterns evolve into a branched network of bands whose edges 
are made up by a series of cusps. These singularities seem to be connected with the 
existence of developable cones and sharp ridges in the free parts of the film, and are not 
yet understood in detail. 
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A A justification of the Donnell equations by formal asymp- 
totic expansion 



The Donnell equations for shells are justified by a formal expansion with respect to a small 
parameter rj proportional to the square-root of the aspect-ratio of the shell. Assuming that 
the transverse displacement is of order tj 2 , we derive these equations from the general non- 
linear equations for elastic shells under finite displacement. The present derivation is 
mainly given for pedagogical purposes as a rigorous proof is available in |31j . Note that 
formal expansions have also been used to justify plate equations in |32j . 



A.l Transformation 

Let us first define the cylindrical basis vectors 





e 1 (6»)= cos# , e 2 (0) = -sin0 . (84) 



This basis is such that §^{9) = e 2 (9), e 2 (#) = —§.i{0), and ^(9) x e 2 {9) = e x - the polar 
axis is the axis X of the shell in reference configuration. 

In the reference configuration, the shell is rolled into a cylinder of radius p. We use 
Lagrangian coordinates (x,y). The position in reference configuration is denoted x, as 
shown in figure [2k, 

— -5- 

Note that the metric associated with the set of coordinates (x, y) is the unit tensor, 
Vx T • = 1. 

The deformed configuration x is defined in terms of the displacement (u x (x, y),u y (x, y),w(x, y)) 
in cylindrical coordinates by equation ([!]), and is also sketched in figure J2ja. 

The strain in the shell is measured using a membrane strain tensor e and a curvature 
strain tensor b. The nonlinear membrane strain e is a 2 x 2 symmetric tensor defined by 
the classical formula: 

F = Vs, £ = Z T -Z, 1=7; (85) 

where the gradient is taken with respect to the Lagrangian coordinates (x,y). The curva- 
ture strain tensor is defined by 

b a p = x >a p(x,y) -N(x,y) (86a) 
where the normal to the shell is defined by 

K(x, y) = -x >x (x, y) x x <y (x, y). (86b) 



The minus sign in the definition of the normal in equation (86b) makes the normal N 
oriented in same direction as the radial vector e 1 (0) of the cylindrical basis. Note that the 
normal is not a unit normal if the current configuration is not developable, e a p ^ 0. Since 



we consider deformations that are almost inextensible, our definition (86a) of the bending 
strain is very close to the geometrically exact definition of the curvature that makes use of 
a unit normal vector, and this introduces a higher-order correction in the thin-shell limit. 



28 



By inserting equation ([!]) into equation (85), one could rederive the fully non-linear 
expression of the 3 independent components e a p in terms of the displacement functions 
(u x ,u y ,w), relevant to the general theory of shells under finite displacements. Similarly, 



by inserting into equation (86), one could derive a fully non-linear expression for the 
curvature strains b a p. 

A. 2 Formal expansion of the membrane and curvature strains 

The Donnell equations are derived from the above set of equations under the assumption 
of a moderate displacement. We consider a formal expansion of the above equations 
with respect to a small parameter tj. The definition of rj given in equation ([9]) will be 
justified. For the moment, it is sufficient to assume ij <C 1. We assume that the transverse 
displacement w scales like p rj 2 , that the in-plane displacement scales like p r/ 3 , and that the 
typical scale for the tangent coordinates x and y is prj. The dimensionless displacement 



u a (x,y) and w(x,y) are defined in equation (11) in terms of the rescaled coordinates 
x = xj {prj) and y = y/(pr)). 

The previous scalings can be justified as follows. Our starting assumption is that the 
deflection w is small, and scales as prj 2 . From this, as w ~ x 2 / p along the curved shell, 
it appears that the natural scale for the tangent coordinates is prj. Balancing the linear 
and non-linear terms in the membrane strain, we have u a , x ~ w 2 x and so the scale for the 
tangential displacement is prj 3 . 

Expanding the membrane strain e introduced in equation (85) with respect to 77, one 
computes 

e a /3{x,y) = 7] 2 e afS (x,y) + C(?? 4 ), (87) 

where the dominant contribution is given by 

- f- -\ fu a ,/3{x,y) + up ia (x,y) \ 1_ 

e a p[x, y) = I V d a y dpy w{x, y) J + - u> iQ (x, y) w^{x, y). (88) 



This result was stated in equation ( 14a) but with the bars omitted. In the right-hand side, 
the derivatives are taken with respect to the rescaled variables, for instance u a B = du x /dy 
when a = x and (3 = y. 



A similar expansion of the curvature strain defined in equation ( 86 ) yields 



77° - 

b aP (x,y) = ^Mx,y) + 0(r7 2 ) (89) 

where the dominant contribution reads 

ba/3 (x, y) = -w ja p (x, y). (90) 

The minus sign in the right-hand side comes from that introduced in the definition of the 
normal N . 

A. 3 Rescaled constitutive equations, elastic energy 

For an isotropic, Hookean (linearly elastic) material, the constitutive laws read, in physical 
units: 

n a p = C {{l-v)e a p + v(tTg)d a p) (91a) 
m a p = D{{l-v)b a p + u{tib)5 aP ) (91b) 
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Figure 10: Equilibrium of a mobile interface G between two domains il 1 and f2 2 . 

where C = Eh/(1 — v 2 ) and D = Eh? /[12(1 — u 2 )] are the stretching and bending moduli, 
respectively. It is convenient to define the slenderness parameter rj by 

1/4 / 1 ^1/2 



n 



D 

Cj 2 



12 p 



(92) 



as we did earlier in equation ([9]). Indeed, this convention makes both the stretching and 
bending moduli C and D effectively equal to one in rescaled units. 
The elastic energy of the shell reads 



(n a/ 3 e a p + m a/3 b a p) dx dy 



and can be rescaled as 



where 



E. 



shell 



p 2 Cr\ E sheih 



E 



shell 



Jla/3 e a p + m a/ 3 b a t 



I dx dy. 



(93) 



(94) 



(95) 



The dimensionless stress can be computed by the dimensionless version of the constitutive 



law, see equation (16). 



As shown in section 2.6, Donnell equations then follow by variational principles from 



the shell energy (95) combined with the definitions (88) and (90) of the membrane and 



curvature strains, and with the rescaled constitutive laws (16). These definitions (88) 



and (90) are often presented as approximation. We have just shown that they follow from 



scaling assumptions and are asymptotically exact for small displacement. 

B Two-dimensional Weierstrass-Erdmann corner conditions 



With the aim to derive the jump conditions at the interface between the adhering part and 
the free part of the shell, we recall the Weierstrass-Erdmann corner conditions in a generic 
two-dimensional setting. We refer to [33] for a detailed presentation. We consider a two- 
dimensional domain $7 that is split in two regions Q± and Q2 meeting along a boundary 
curve G = Sli n SI2, as depicted in figure 10 The unknowns are the functions £ a (a;, y) 



where a is an index, such as the component of the displacement in the case of an adhering 
shell. Each region S7j, i = 1, 2 is associated with a specific Lagrangian (also called energy 
functional), 



fii(G) 



£ l (£,a, dx dy. 



(96) 



30 



The boundary may evolve and we use the boundary curve G as an unknown: the domains 
fij are reconstructed in terms of the curve G, hence the notation f2,(G). 

We are interested in the conditions that make the total energy E — E\ -j- E2 minimum, 
and in particular in the conditions associated with the motion of the interface G. The 
variation 5E of the total energy is made up of a term arising from the variation of the 
functions £, a (x, y) inside each domain, and from a term associated with the motion of the 
free boundary, 

^Amg) V<96* dU,p J Jg 

Here the partial derivatives denote functional derivative, d£ is the element of length along 
the boundary G, and 8Gd£ is the signed area swept by the free boundary, counted posi- 



tively when the region Q2 grows while the region Q\ shrinks, as shown in figure 10 

We consider the contribution 5Eq to SE that collects all terms written as integrals 
along the boundary curve G. These terms yield the jump conditions associated with 
the equilibrium of the boundary, while the other terms contribute to the Euler-Lagrange 
conditions of equilibrium inside each subdomain. This 5Eq is made up of the last term in 



equation (97), and of a term coming from the integration by part of the term proportional 

to 5^ a)l 3, 

5E G = j G (V] 5G + 5U Np} d£. (98) 

Here the double bracket denotes the jump, [/*] = f 2 — f 1 and N = (N x ,N y ) is the vector 
normal to the boundary G, directed towards region 1 as in the figure. Along the common 
boundary G, this N_ is equal to the outward normal N 2 with respect to the domain Q2, 
and is opposite to the outward normal iV 1 to domain Q\. 

We shall assume that the unknown £ a is prescribed to be continuous across the bound- 
ary, as happens with the components of the displacement and with the slope q a = w a of 
an elastic shell, 

KJ = 0. (99) 

This continuity relation can be differentiated in a frame moving along with the boundary. 
This yields, for an arbitrary perturbation of the boundary and of the function, 

lSU + [^p}NpSG = (100) 

First, consider perturbations leaving the boundary unchanged, 5G = 0. Then equa- 



tion ( 100 ) shows that 5£ a is continuous across the boundary as [<5£q,] = 0. In the expression 



for 5 Eg given in equation (98), the first term cancels when the free boundary is at rest, 

operator, yielding the condition 

1^6^ = 0. (101) 



and 5£ a can be factored out of the jump operator, yielding the condition 



In the particular case of a function £ 7 whose values are prescribed in the domain fii, 
we have <5£* = 0. The continuity condition (101) then yields + [£ 7l /3'] Npi 8G = 0. 



Cancelling the variation in equation ( 98 ) then yields 



dC 2 

if <5£ 7 is prescribed on fix, then C 2 - C 1 - — — [£ 7 «/] Np> Np = 0. (102) 
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